Scaling relations in anisotropic superconductors with strong pair-breaking 
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Following the work of Abrikosov and Gor'kov on the pair-breaking effects, one can derive the 
temperature dependences of the electronic specific heat C s /T = 7' + 1-iT 2 (with the jump at the 
superconducting transition AC oc for materials with zero Fermi surface average of the order 
parameter (A) = (e.g. d-wave) or for those with (A) <C A maa; (e.g., ±s of iron-pnictides) in 
the presence of strong pair-breaking. Moreover, the London penetration depth satisfies A~ 2 = 
A^ 2 (l — T 2 /T 2 ) (or A — Ao = /3T 2 at low temperatures) and the slope of the upper critical field 
near T c is H' c2 oc T c . Remarkably simple relations between these at first sight unrelated quantities 
take place: fJ-X^T^ /\H' c2 \ = 3<f>o/8n 2 and AC /3 2 T£/\H' c2 c \ = </>o/167r 2 are universal constants. The 
prediction is checked on two samples of Ba(Fei_ a; Co :E )2As2 and on CeCoIns for which the data 
needed are available. 

PACS numbers: 74.20.-z, 74.20.Rp 



A remarkable scaling relation between the specific heat 
jump in iron-pnictides and their critical temperatures 
has recently been established: AC oc T c 3 The same 
relation, in fact, follows from the theory developed by 
Abrikosov and Gor'kov (AG) for isotropic superconduc- 
tors with magnetic impurities |2j of a concentration close 
to the critical for which T c <C T c0 , the critical tempera- 
ture in the absence of spin-flip scattering. This prompted 
the author to derive the Ginzburg-Landau (GL) equa- 
tions for the case of a strongly anisotropic order parame- 
ter with a strong pair-breaking [3| , that made it possible 
to obtain the jump 



AC — C s C n 



l§ir 4 k%N{0)T 2 + 
3?i 2 (3(tt 4 ) - 2) 



(1) 



Here, A^(0) is the total density of states at the Fermi level 
per one spin. The order parameter was taken in the form 
A(r, T;kp) — *(r, T) fl(k F ) where fl(k F ) describes the 
variation of A along the Fermi surface and is normalized 
to have the average over the whole surface (f7 2 ) = 1. 

Another assumption employed in the derivation of 
Eq. fTJ was that (A) = (£1) = (which is true for the 
d-wave) or that (A) <C A max (presumably true for iron- 
pnictides with ±s order parameter). 

The scattering in the Born approximation was charac- 
terized by two scattering times, the transport r respon- 
sible for the normal conductivity and T m for processes 
breaking the time reversal symmetry (e.g., spin- flip): 



1/t+ = 1/t + 1/t„ 



(2) 



This is of course a gross simplification. For multi-band 
Fermi surfaces one may need more parameters for vari- 
ous intra- and inter-band processes; these, however, are 
hardly controllable and their number is too large for a 
useful theory. 

It is well-known that the formal AG scheme for uncor- 
rected magnetic impurities |2| applies to various situa- 
tions with different pair-breaking causes, not necessarily 



the AG spin- flip scattering |4( . A particular example rele- 
vant for iron-pnictides is antiferromagnetic superconduc- 
tors in which even non-magnetic impurities cause pair- 
breaking, see e.g. review 

Another quantity obtained in [3] from the GL equa- 
tions is the slope of H C 2 at T c : 



dH, 



c2,c 



dT 



'37i 2 (fi3«3) 



(3) 



for the field along the c axis of a uniaxial crystal (c/>o is 
the flux quantum). The scaling relation H' c2 oc T c holds 
across the whole family of iron-pnictides, albeit with a 
considerable scatter; a number of data points are col- 
lected in Q to support this statement. There are rea- 
sons for this scatter: different criteria in extracting H c2 
from resistivity data, unavoidable uncertainties rooted in 
sample inhomogeneities in determination of T c and the 
slopes of H C 2(T) near T c , possible differences in Fermi ve- 
locities and the order parameter anisotropics, to name a 
few. Since the model employing only two scattering pa- 
rameters for multi-band iron-pnictides is a far-reaching 
simplification, one can expect the model to work quali- 
tatively at best. 

The third piece of evidence in favor of the strong pair- 
breaking in iron-pnictides is the low temperature behav- 
ior of the London penetration depth X(T). Evaluation of 
A(T;r, r m ) for arbitrary r's and arbitrary anisotropy of 
A is difficult analytically. However, for a strong T c sup- 
pression, as was shown by AG Q , the formalism used for 
derivation of the GL equations near T c applies at all tem- 
peratures. Physically, this is because the pair-breaking 
suppresses the order parameter so that the expansion in 
powers of A and its derivatives can be done at any T. 
The calculation then proceeds in a manner similar to 
that near T c Within the same formal scheme that 
gave Eq. (UJ for the specific heat jump and Eq. Q for 
the H C 2 slopes, one obtains for the order parameter in 



2 



the field-free state: 



* 2 



2n*k%(T?-T*) 
3(ft 4 ) - 2 



(4) 



this reduces to the AG form for O = 1. The free en- 
ergy difference between the normal and superconducting 
states is given by Q: 



F — F 



3r(3(0 4 ) -2) 



(5) 



Hence, one can evaluate not only the jump AC at T c , but 
the full temperature dependence C S {T) = -Td 2 F s /dT 2 : 



£ = ( 7 - M T c 2 /3)+/xT 2 . 



8^ 4 fc 4 A(0)r 2 
h 2 (3(n 4 } -2) 



(6) 
(7) 



Here, jT is the electronic heat capacity of the normal 
phase. The dependence of the type ^ has recently been 
reported in Refs. 

for Ba(Fei_ 3; Co 3; )2As2. 
Within the same model one can also calculate the Lon- 
don penetration depth Q: 



16n 3 e 2 N(0)k 2 B T 2 
c 2 h 2 (3(n 4 ) -2) 
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The low T behavior of AA ab 
uniaxial material is given by: 



Kb(T) - A ab (0) for a 



AA. 
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The power-law behavior of AA ab oc T n with n ~ 2 has 
been seen in many iron-pnictides Moreover, the pre- 
factor rj/T 3 extracted from the data varies from one ma- 
terial to another approximately as 1 /T c 3 that further sup- 
ports the pair-breaking scenario. One readily obtains for 
T = 0: 



Aah(O) = 2 V /T C 



(10) 



The relations ©,(0) for the specific heat, flSJ for 
the slope of H c2 , and (|8|)-([10"]) for the penetration 
depth contain material parameters: the density of states 
N(0), Fermi velocities Vi, the parameter f2 of the gap 
anisotropy, and the combined scattering time r+. One 
can express the ratio (3(£1 4 ) — 2)/rl in terms of the co- 
efficient /i of Eq. (|6|) and substitute the result in the ex- 
pression for A(0), thus establishing a relation between the 
behavior of the specific heat and the penetration depth. 
(v 2 n 2 ) can in turn be expressed in terms of the H c2 slope 
of Eq. ([3]) . Surprisingly, this simple algebra yields 



MAg b (0)T c 3 _ 30o 
8tt 2 



(11) 



with a universal quantity on the right-hand side (RHS). 
Four quantities on the LHS are measured in indepen- 
dent experiments so that this result might be used as a 
stringent test of the pair-breaking scenario. It is worth 
reiterating that this relation is obtained for nearly zero 
(A)'s in materials with a strong pair-breaking. 

To check the scaling relation (fTTj) one can turn to data 
on Ba(Fei_ 2; Co :l; )2As2 with x = 0.075 and T c = 21.4 K 
for which the behavior of the electronic specific heat has 
recently been reported 0] • The phonon background was 
evaluated by measuring the specific heat of an overdoped 
material with x — 0.153 for which superconductivity is 
suppressed. This procedure yields C s (T) of the form ([6]) 
with 



H « 0.12mJ/molK 4 w 20erg/cm 3 K 4 



(12) 



There is still a considerable uncertainty in experimen- 
tal values for zero-T penetration depth which in var- 
ious studies with different techniques is found to be 
A ab (0) « (270 ± 50) nm 0, [l(J. Taking H' c2 c « 2.5 T/K, 
see [3JI, one estimates: 



\H' c2 J 



(3.7-7.8) x 10" 9 Gcm 2 



(13) 



This is of the same order of magnitude as the RHS of 
Eq. (JTIJ: 3(/>o/87r 2 w 7.6 x 10- 9 Gcm 2 . 

Similar specific heat results have recently been re- 
ported for the same compound with x — 0.08, albeit 
with a different way to subtract the lattice contribution 
ll|. The low temperature behavior of C s /T is linear 
in T 2 with approximately the same value of the slope 
d(C s /T)/d(T 2 ) as reported in 0. 

Since other iron-pnictides comply with the scaling rela- 
tions ©, ©, and J3J (see 0, Q and 0), they are good 
candidates to obey the relation (fTTj) . Unfortunately, for 
a precious few of these compounds measurements of all 
quantities involved in Eq. pip were done on the same 
samples. A good example is Ref. [|| where C s /T of a few 
crystals of Ba(Fei_ K Co :c )2As2 was found to be linear in 
T 2 (C (T) in the field of 9 T was taken as the phonon con- 
tribution and was subtracted). However, neither H' c2 (T c ) 
nor A(0) are available for these samples. In this situation 
one can use a crude estimate of slopes H' c2 (T c ) based on 
comparison of the scaling relation ([3]) with data on many 
iron-pnictides 0, H' c2 (T c ) w (0.2 T/K 2 ) T c to obtain val- 
ues listed in Table 1. Then, with the help of Eq. (fTTj) one 
can calculate the penetration depth. The values of A(0) 
so obtained and given in the Table are quite reasonable. 

As mentioned, the results discussed above hold for 
strong pair-breaking in materials with the order param- 
eter averaged over the Fermi surface close to zero. This 
is exactly the case for the d-wave compounds. Hence, 
the underdoped cuprates with a strongly suppressed T c 
may also be among materials satisfying the above scaling 
formulas and the scaling relation (|lll) . in particular. It 
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TABLE I: The penetration depth A(0) calculated with the 
help of Eq. (fTTjl for Ba(Fei_^Co x )2 As2. The slopes /i — 
d(C s /T)/d(T 2 ) are extracted from Fig. 3b of Ref ||. T c is 
taken as corresponding to zero resistance. H' c2 (T c ) is evalu- 
ated as explained in the text. 



immediately follows: 



x 


T c K 


/imJ /K 4 


H' c2 (T c ) T/K 


A(0)nm 


0.06 


14.5 


0.8 


2.9 


240 


0.07 


19.2 


0.67 


3.8 


200 


0.08 


23.9 


0.73 


4.8 


150 


0.12 


17.3 


0.75 


3.4 


200 


0.15 


7.0 


0.77 


1.4 


500 



is worth noting that the scaling (fTUf is supported by the 
surface resistance [l3| and optical data 14 1 for a series of 
samplesof YBa2Cu306+£c with T c varying from 3 to 17 K. 

Another example where the model developed here may 
work is CeCoIns. This is a clean heavy-fermion d-wave 

2.3 K; in the normal 



superconductor 15|, [16J with T 
phase, the material is a paramagnet [13] • Note that quan- 
tities entering the relation (fTTj) are either for zero-field (/i 
and A(0)) or for temperatures near T c (H' c2c ), so that the 
domain of interest here is not affected by complications 
related to paramagnetic constrains or to possible Fulde- 
Ferrell-Larkin-Ovchinnikov phase. 

The temperature dependence of the penetration 
depth of CeCoIn 5 [l2| is well described by A = 
A(0)/a/I -T 2 /T c 2 (with A « 360 nm) [T3] as is should 
be for a strong pair-breaking, Eq. ((5]). Moreover, C s /T 
(where C s is the electronic part of the specific heat) 
is close to being linear in T 2 at low temperatures [lSj] 
in agreement with Eq. (|6]). Reading the data on C s /T 
on Fig. 2b of [Hj], one can estimate the coefficient /i of 
the T 2 part as \i sa 3.5 x 10 4 erg/cm 3 K 3 . The slope of 
the upper critical field is known too: H' c2 c ps 11.5 T/K. 
Then, the LHS of the scaling relation (ITTj) is estimated 
as 4.8 x 10 _9 Gcm 2 , reasonably close to the universal 
number 7.6 x 1CP 9 G cm 2 on the RHS of this relation. 

The scaling relation involves two quantities, the 
slope n = d(C s /T)/d(T 2 ) and the zero-T penetration 
depth A(0), both of which are difficult to access exper- 
imentally (the first because of the necessity to subtract 
the phonon contribution from the measured C/T, the 
second because usually only the deviation of A(T) from 
A(0) is measured). Fortunately, one can avoid these diffi- 
culties using another scaling relation. Writing Eq. © as 
AA a6 = j3T 2 with (3 = 7?/T c 3 one readily finds: 

AC/3 2 T 4 6 a 



\H' c2 J 



16tt 2 



1.27 x lCT a Gcrr/ 



(14) 



Measurement of the jump AC does not require the 
phonon contribution be subtracted and determination of 
j3 = d(AX) / d(T 2 ) does not require knowledge of A(0). 

It is worth noting a couple of other scaling relations 
which might be useful. From expressions ([T]) and ([7]) 



2T 2 



Similarly, Eqs. ^ and (JTO]) yield 
2T 2 



A(0) = ^f- AA(T). 



(15) 



(16) 



A word of caution here is in order. The scaling re- 
lations discussed are obtained for superconductors with 
(A) « and in the limit of strong pair-breaking, strictly 
speaking for T c — > 0. It is shown that in this limit C s /T 
is linear in T 2 and A A is proportional to T 2 . For sys- 
tems discussed above, this is only approximately so. E.g., 
for CeCoIn5 the actual fit of the experimental C s /T [18j 
to the form /ii + pT v yields v k, 2.3 indicating that 
the material is not quite in the limit of a strong pair- 
breaking. Similarly, the data @ yield the power-law for 
AA a fc oc T n with n close but not exactly equal to 2, so 
that the iron-pnictides examined are also not exactly in 
the strong pair-breaking limit. This, along with common 
experimental uncertainties, is, probably, the main source 
of deviations from the scaling relations proposed. 

To summarize, scaling relations are found involving the 
specific heat, the penetration depth, and the upper criti- 
cal field for superconductors with a strong pair-breaking 
and the averaged over the Fermi surface order parameter 
close to zero (as is the case for d-wave or for ±s symme- 
tries with impurities). The relation are checked on a few 
iron-pnictides and on CeCoIn 5 . 
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